Characterizing top gated bilayer graphene interaction with its environment by Raman 

spectroscopy 
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In this work we study the behavior of the optical phonon modes in bilayer graphene devices by 
applying top gate voltage, using Raman scattering. We observe the splitting of the Raman G band as 
we tune the Fermi level of the sample, which is explained in terms of mixing of the Raman {Eg) and 
infrared {Eu) phonon modes, due to different doping in the two layers. We theoretically analyze our 
data in terms of the bilayer graphene phonon self-energy which includes non-homogeneous charge 
carrier doping between the graphene layers. We show that the comparison between the experiment 
and theoretical model not only gives information about the total charge concentration in the bilayer 
graphene device, but also allows to separately quantify the amount of unintentional charge coming 
from the top and the bottom of the system, and therefore to characterize the interaction of bilayer 
graphene with its surrounding environment. 

PACS numbers: 02.20.-a, 78.30.-j, 78.67.-n 



Bilayer graphene has attracted a lot of attention re- 
cently because of its special low energy electronic dis- 
persion, in which a tunable band gap can be opened by 
application of a transverse electric field [H-ll]. Such de- 
vice is desirable for low energy photo-emitters and de- 
tectors possessing a high tunability by the control of 
charge concentrations on the graphene layers. Recent 
experimental demonstration of this tunable band gap in 
bilayer graphene was based on the absorption measure- 
ments in the infrared region or by electric transport 
measurements [3, S]- However the tunable band gap bi- 
layer graphene device operation can be greatly influenced 
by the surrounding environment. Typically, uninten- 
tional doping charges coming from the top and the bot- 
tom of the system can accumulate on bilayer graphene, 
giving rise to an unintentional electric field which deter- 
mines a non-homogeneous doping between the layers and 
the opening of a band gap in the band structure, without 
any applied electric field @ . In this work we use Raman 
spectroscopy to monitor the unintentional charge coming 
from the top and the bottom of the system, which gives 
information on the electrostatic environment of the sam- 
ple and which helps to characterize the bilayer devices 
for further applications. 

The band gap opening and tunability in bilayer 
graphene is based on the application of an electric field 
E perpendicular to the layers, given hy E = {ritop — 
nbot)\e\ I {'^^o) where ritop and ribot are the charge car- 
riers coming from the top and the bottom of bilayer 
graphene, respectively. Raman spectroscopy has already 
shown to be a fast and non-destructive tool to character- 
ize graphene samples [13, [HI and doping effects [l^ - [l6| , 
however no carefully analysis has been done to demon- 
strate the effect of non-homogeneous doping in bilayer 
graphene devices. Recent theoretical calculations made 
by Gavaet al. @ suggest that from the analysis of the 



Raman spectra of gated bilayer graphene it is possible 
to quantitatively identify the amount of non-intentional 
charges coming from the atmosphere and from the sub- 
strate and to characterize the electrostatic environment 
of few-layers graphene. In this work we study the de- 
pendence of the G band of bilayer graphene on the gate 
voltage. From the direct comparison between the exper- 
imental and the theoretically simulated Raman spectra, 
and from the analysis of the positions, full width at half 
maximum (FWHM) and relative intensities of the two 
Raman peaks as a function of the electron concentra- 
tion, we were able to estimate the charge unintentionally 
accumulated on the device from the environment. 

Fig. [Ija) shows the bilayer-graphene field-effect tran- 
sistor (FET) used in the experiment. Graphene samples 
were produced by micro-mechanical cleavage of graphite 
and deposited on Si covered with 300nm of Si02- Top 
gating was achieved by using a polymer electrolyte con- 
sisting of polyethylene glycol (PEG) and NaC104 with 
ratio concentration of 1:0.25, and the gate voltage was 
applied between a gold electrode in contact with the 
graphene layer and a platinum wire electrode inserted 
in the electrolyte (see the schematic setup in Fig. [TJb)). 
The contacts were made by optical lithography. The Ra- 
man measurements were done in the back scattering con- 
figuration at room temperature using a triple monochro- 
mator spectrometer (DILOR XY) using 2.41 eV as exci- 
tation laser energy. The spot size of the laser was '^1 ^m 
using a 80 X objective and the laser power was kept at 
1.4 mW. 

The interface between graphene and polymer elec- 
trolyte has been shown experimentally to behave like 
a double layer capacitor of thickness in the order of 
nanometers . Therefore, the geometric capacitance of 
the electrolyte is very high compared to bottom gate de- 
vices where the thickness of the dielectric is much larger 



2 




Pt-wire 



FIG. 1: (a) Optical microscope image of the graphene sample 
before application of the polymer electrolyte, (b) Schematic 
illustration of the device and the experimental setup. 



(typically 300 nm). 

Capacitance measurements of the polymeric elec- 
trolyte used in the experiment were performed by 
Impedance Spectroscopy with frequency analyzer AU- 
TOLAB PGSTAT30 by using a symmetrical cell with 
two Au electrodes and a polymer electrolyte layer. The 
impedance measurements for the electrolytes were car- 
ried out with frequency ranging from 50 kHz to 0.5 Hz 
at V with 10 mV amplitude, and we obtained the value 
of 1.5 X 10"^ F cm~^. However, the shape and the thick- 
ness of the electrical double layer depends on the specific 
surface at contact with the electrolyte. Therefore, the 
measured value of Cg only gives the order of magnitude 
of the geometric capacitance of the electrolyte in contact 
with bilayer graphene. 

In the case of bilayer graphene with Bcrnal AB layer 
stacking, both the electronic and phonon bands split into 
two components with special symmetries [17[. The 
phonon mode of monolayer graphene splits into two com- 
ponents, associated with the symmetric (S) and anti- 
symmetric (AS) displacements of the atoms in the two 
layers with respect to inversion symmetry fisj . The S 
and AS modes belong to the two double degenerated rep- 
resentation Eg and E„, respectively [ol. [l8H20i |. The E„ 



mode is not Raman active and, therefore, the G band of 
isolated bilayer graphene is composed of only one peak. 
However, when the two layers of bilayer graphene have 
different charge carrier concentration, induced by the ap- 
plication of an external gate voltage, the inversion sym- 
metry of bilayer graphene is broken, lowering the sym- 
metry of the system. As a consequence of the induced 
asymmetry between the two layers, the two S and AS 
modes are mixed, the two new eigen-modes have the Ra- 
man active S component, and therefore two pea ks are 
observed in the G band of bilayer graphene (ol. Il8l - [20j . 

In Fig. HJa) (red dashed curves) we show the experi- 
mental Raman spectra taken with the application of top 
gate voltage (Vg) from -1.50 to 1.00 V. The G band split- 
ting into two components Gh and G; (higher and lower 
frequency peak, respectively) can be clearly observed for 
Vg below -0.6 V. In order to compare the experimen- 
tal spectra and theoretical calculations, we converted Vg 
into n using the expression /3Vg = (n — no)e where the 
total capacitance /3, which includes the quantum capac- 
itance (Cg) and geometrical capacitance (Cg) [13], and 
the intrinsic doping at the zero gate uq are used as fit- 



ting parameters. Moreover, by the comparison between 
experimental and theoretical results we can estimate the 
charges unintentionally adsorbed, at zero gate, from the 
top and bottom of the device, n^^p and n^^^ respectively. 



These quantities are related to ng by rto 
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and 



therefore we only used n^^^ as additional fitting parame- 
ter. Finally, the theoretical FWHM F*'* calculated as a 
function of n and n^^^ is given by electron-phonon and 
an-harmonic phonon- phonon interaction Therefore, 
in order to take into account other factors determining a 
finite lifetime and neglected in the calculations, we used 
in the fitting procedure a parameter F°, independent on 
the total charge n and equal for the two peaks, related 
to the total FWHM by F = F*'' + F°. 

In Fig[2Ka) we show the comparison between the exper- 
imental spectra (red dashed curve) and theoretical one 
(black continuous curve), obtained using the parameters 
discussed below, for different Vg. The fit is performed 
computing the square of the difference between the ex- 
perimental and theoretical spectra, averaged over all the 
measured Raman range and over different Vg. We con- 
sidered Vg in the range of ±0.5 V. This choice is moti- 
vated by the fact that for large values of Vg the linear 
relation between gate voltage and charge could be mod- 
ified, and charges from the electrolyte could accumulate 
on the bottom of the sample, making the fit results less 
reliable. The theoretical spectra is obtained as the sum of 
two Lorentzians, for the G^ and G/ frequency peaks. The 



two Lorentzians are centered in 



with FWHMs F 



h/l 



Fl"/, -I- F°, and with area Ih/i, where w/j/;, F*^'^j, and 



^h/; -'^ '^^^'^ wiLU area i/^/i 
arc computed as a function of n and nj^^^j. In par- 
ticular, we used two different parameters, f3~^ and /3~, 
for positive and negative Vg, which induce positive and 
negative n (i.e. electron and hole doping charge), re- 
spectively. The values for the parameters used in the fit 
(i.e. /S"*", /3~, no, n^^^, and F") are varied in uniform and 
dense grids. The best fit is obtained for /S"*" = 3.7 and 



4.6 X 10"'' F cm- 



no 



-0.15 X 10^3 c^-^ 



^hot 



0.0, and F° = 7.5 cm~^. Notice that the agreement 
between the experimental and simulated spectra is excel- 
lent in the range -0.5 to 0.5 V. The slight shifts out of this 
range can be ascribed to a possible hysteresis in the ex- 
periment, and to the fact that we did not consider in our 
model the expected jump of the quantum capacitance 
(Cq) when, increasing (decreasing) the Fermi level, we 
reach the second conduction band (penultimate valence 
band) in bilayer graphene (l3 |. The different values of 
Z?"*" and /3~ can be ascribed to the different mobilities of 
the positive {Na~^) and negative (ClO^) ions. 

All the G band experimental spectra of Fig.[2l^a) were 
fitted using two Lorentzians in order to extract the fre- 
quencies and full width at half maximum (FWHM) of 
two components Gh and G/, as well as the relative Ra- 
man intensity, i.e., the ratio between the Raman intensi- 
ties of the modes with highest frequency (1^) and lowest 
frequency (I/). Figs.[HJb) and (c) shows, respectively, 
the dependence of the Gh (red dots) and G/ (black dots) 
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FIG. 2: (a) Comparison between experimental (red dasiied 
line) and theoretical (black continuous line) spectra for dif- 
ferent Vg as obtained from the direct fit of spectra; (b-c-d) 
Frequency uj^/i, FWHM F^y; and ratio of I/i/; as obtained fit- 
ting the experimental spectra with two Lorentzians and com- 
pared with theory. The dots are the experimental data and 
the full curves are the theoretical calculations. 

frequencies and FWHM F as a function of the electron 
concentration. The full lines are the theoretical calcula- 
tions. The dependence of the frequency of the Gh and 
G; Raman peaks (Fig. ^b)) is well described by the 
calculation of the phonon self energy as a function of 
charge concentration. The distinct behaviors of the Gh 



and G/ is qualitatively explained by the distinct electron- 
phonon couplings of these modes. While Gh blue-shifts 
with charge carrier concentration due to inter-band tran- 
sitions, the G; mode redshifts due to intra-band transi- 
tions when the Ej? is changed. For the FWHM depen- 
dence, while F; does not change with n, the Th is max- 
imum near n = and minimum for values of n corre- 
sponding to values oi Ep larger than half of the phonon 
energy, as has been observed before in both mono and 
bilayer graphene It's worth to mention that the 

scattered data points for the FHWM are mostly caused 
by charge carrier fluctuation during the measurement, 
where a hysteresis of the charge neutrality point is found 
by sweeping the gate voltage up and down. In Fig. [2][d) 
we plot the dependence of the Ih/h as a function of n. 
The quantity Ih/h shows a minimum value between n^ 
-1 to 1 X e 10^'^ cm~^, and increases more strongly for 
positive carrier concentration. 

In summary, a detailed analysis of the G band of 
top gated bilayer graphene is presented. We observed 
that, unlike in the unbiased case where the G Raman 
band is composed by only one peak, the gate voltage 
breaks the inversion symmetry and the G band splits 
in two modes, that are combinations of the symmet- 
ric and anti-symmetric modes of the unbiased bilayer 
graphene. We analyze the dependence of the frequency 
and the relative intensities of the peaks with higher and 
lower frequency as a function of the electron concentra- 
tion and we compared the experimental results with the- 
oretical calculations. From this comparison, we could es- 
timate the unintentional carrier concentration adsorbed 
on the device, at zero gate, from the substrate, and 
from the electrolyte, n^^p, and we found n^^^ = 0.0 and 
rijQp ~ —0.15 X 10"'^'^ cm~^, showing that Raman spec- 
troscopy is a powerful technique to study the electrostatic 
environment of graphene. 
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